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ABSTRACT 

We present 5 — 38/im mid-infrared spectra at a spectral resolution of R « 65 — 130 of a large sample 
of 22 starburst nuclei taken with the Infrared Spectrograph IRS on board the Spitzer Space Telescope. 
The spectra show a vast range in starburst SEDs. The silicate absorption ranges from essentially no 
absorption to heavily obscured systems with an optical depth of rg.8 Mm ~ 5. The spectral slopes can 
be used to discriminate between starburst and AGN powered sources. The monochromatic continuum 
fluxes at 15/xm and 30//m enable a remarkably accurate estimate of the total infrared luminosity of 
the starburst. We find that the PAH equivalent width is independent of the total starburst luminosity 
Lir as both continuum and PAH feature scale proportionally. However, the luminosity of the 6.2/im 
feature scales with Lm and can be used to approximate the total infrared luminosity of the starburst. 
Although our starburst sample covers about a factor often difference in the [Nelll] / [Nell] ratio, we 
found no systematic correlation between the radiation field hardness and the PAH equivalent width 
or the 7.7/mi / 11.3/mi PAH ratio. These results are based on spatially integrated diagnostics over 
an entire starburst region, and local variations may be "averaged out" . It is presumably due to this 
effect that unresolved starburst nuclei with significantly different global properties appear spectrally 
as rather similar members of one class of objects. 

Subject headings: Telescopes: Spitzer, Galaxies: ISM, starburst, star clusters ISM: dust, extinction, 
HII regions 
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1. INTRODUCTION 

Many nearby galaxies show dramatically increased 
rates of star formation compared with the Milky Way. In 
such "Starburst" galaxies (e.g., Weedman et al. (1981)), 
the primary energy source is driven by high nuclear star 
formation rates, rapid nuclear gas depletion time-scales 
and high supernovae rates. Such starburst systems often, 
though not exclusively, occur in interacting and colliding 
systems. Since collisions and interactions are believed 
to be a fundamental part of the evolution of galaxies 
by hierarchical growth, the full characterization of star- 
burst galaxies is of the great importance in measuring 
and quantifying the global history of star formation over 
cosmic time. Determining the average mid-infrared (mid- 
IR) spectral properties, and the range of observed behav- 
ior within the starburst class at low redshift is vital for 
interpreting spectra of higher redshift IR sources, pro- 
viding complementary spectral templates to those paral- 
lel Spitzer studies of the more extreme ULIRG systems 
(e.g., Armus et al. (2006b)) and AGN (e.g., Weedman et 
al. (2005)). 

The term "starburst galaxy" is commonly used to de- 
scribe an apparently well-defined class of objects, al- 
though starbursts can be found in the most diverse condi- 
tions, ranging from low-pressure dwarf galaxies to high- 
pressure nuclear starbursts. Their observed properties 
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are expected to depend on numerous parameters such as 
the initial stellar mass function (IMF), the duration and 
epoch of the individual starburst (s), the metallicity of 
the ISM, the size and distribution of the dust grains, the 
strength of the magnetic fields, gas pressure and temper- 
ature of the ISM, galactic shear, total luminosity, and 
total mass. Furthermore, nearby starbursts, for which 
high resolution imaging is possible, have revealed com- 
plex substructures - in both stellar distributions and ISM 
- ranging from ultra-compact Hn regions (UCHIIR) to 
large complexes of super star clusters (SSC), suggesting 
small-scale variations of the observables across a star- 
burst region. 

We use the low resolution mode of the Infrared Spec- 
trograph 3 (IRS) (Houck et al. 2004) on board the Spitzer 
Space Telescope (Werner et al. 2004) to observe the cen- 
tral regions of 22 starburst galaxies. Our objects repre- 
sent a sample of "classical" starbursts for which a wealth 
of literature exists. The sample includes both purely 
starburst and starbursts with weak AGN activity (as 
determined from X-ray, optical, or radio observations). 
The summary in Table 1 lists the observed targets, their 
general properties, the classifications we adopt, and the 
references from which they are derived. The continuous 
5 — 38/im IRS spectra include the silicate bands around 
10/im and 18/im, a large number of PAH emission fea- 
tures, and information on the slope of the spectral con- 

3 The IRS was a collaborative venture between Cornell Univer- 
sity and Ball Aerospace Corporation funded by NASA through the 
Jet Propulsion Laboratory and the Ames Research Center. 
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tinuum. 

Numerous mid-IR studies of starbursts have been con- 
ducted with ISO-SWS, ISOCAM, or ISOPHOT, see for 
instance, Rigopoulou et al. (1996), Lutz et al. (1998), 
Rigopoulou et al. (1999), Dale et al. (2000), Helou et al. 
(2000), Laurent et al. (2000), Sturm et al. (2000), Thorn- 
ley et al. (2000), Charmandaris et al. (2001), Forster 
Schreiber et al. (2003), Lu et al. (2003), Verma et al. 
(2003), Tacconi-Garman et al. (2005), and Madden et 
al. (2006). An overview of many ISO results is given in 
Genzel & Cesarsky (2000). While the ISO observations 
provided great new insights in the spectral characteris- 
tics of individual starbursts, the sample of continuum 
spectra remained rather small or was limited to shorter 
wavelengths or narrow bandwidth scans of the strongest 
emission lines. 

In this paper we will address the question of the mid- 
IR homogenity of the classical starburst class 4 , and at- 
tempt to investigate how specific spectral features (espe- 
cially the mid-IR PAH bands) vary with the total UV 
continuum flux, UV hardness ratio and dust extinction 
within the starburst nucleus. We will also investigate 
how the shape of the continuum depends on the lumi- 
nosity source, and if the total luminosity affects the ob- 
served spectral shapes, i.e., to what extent starbursts 
can be scaled up. We investigate the role of dust, and 
how well PAH emission correlates with the rate of star 
formation. 

The outline of the paper is as follows: First we give a 
detailed description of the observations and the data re- 
duction and calibration. In Section 3 we discuss how the 
relevant spectral features (SED, PAHs, silicate features) 
have been measured. The main focus will be on the dis- 
cussion of the numerous results in section 4, followed by 
a summary. We note that our sample has also been ob- 
served with the IRS high resolution modules, revealing 
the large, comprehensive zoo of strong and faint fine- 
structure lines in the 10 — 38/im wavelength range. This 
analysis is complementary to our above science goals and 
will be presented in a subsequent paper by Devost et al. 
(2006). 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations 

We observed all targets with the two low-resolution 
modules (R « 65 - 130) of the IRS. The slit widths are 
about 3.6" from 5 — 15/im and 10.5" from 15 — 38/im. 
The observations were made within the first year of the 
Spitzer Space Telescope mission (see Table 2) as part of 
the IRS guaranteed time observing program. The data 
were taken using standard IRS "Staring Mode" Astro- 
nomical Observing Templates (AOT). In most cases, a 
high accuracy IRS blue peak-up, offsetting from a nearby 
2MASS star, was performed to achieve the intended 
pointing accuracy. The central coordinates of the slits 
for these observations were derived from 2MASS images. 
Table 2 lists the observing parameters for all targets. 
Fig. 1 shows the slit positions relative to the galaxies. 

2.2. Data Reduction 

4 In this paper the terms 'starburst', 'starburst galaxy', and 
'starburst nuclei' all refer to the central, sub-kiloparsec regions of 
galaxies with significantly enhanced starburst activity 



The data were pre-processed by the Spitzer Science 
Center (SSC) data reduction pipeline version 11.0 (SOM 
2004) (except for NGC3256, which was processed with 
version 12). To avoid uncertainties introduced by the 
flat-fielding in earlier versions of the automated pipeline 
processing we started from the two-dimensional, unflat- 
fielded data products, which only lack stray light cor- 
rection and flat fielding. These products are part of the 
"basic calibrated data" (BCD) package provided by the 
SSC. The various steps of the data reduction followed a 
recipe that has been developed by the IRS Disks team 
and tested on a large amount of Galactic and extragalac- 
tic spectra. 

We used the Spectral Modeling, Analysis, and Reduc- 
tion Tool (SMART) version 5.5.1., developed by the IRS 
team (Higdon et al. 2004) to reduce and extract the spec- 
tra. First, we median-combined the images of the same 
order, same module and same nod position. Then we dif- 
ferenced the two apertures to subtract the "sky" back- 
ground, which is mainly from zodiacal light emission. 
We have done so by subtracting the spatially offset first 
and second order slits from each other: SL1 from SL2 
and vice versa for the SL (short wave, low resolution) 
module, and LL1 from LL2 and vice versa for the LL 
(long wave, low resolution) module. We note that this 
approach assumes that the emission from the target is 
not farther extended than the angular distance between 
the two corresponding subslits of 79" in SL. 

We extracted the spectra using a column width which 
increases - like the instrumental point-spread function 
- linearly with wavelength. The extraction width is set 
to four pixels at the central wavelength of each subslit. 
The spectra were flat-fielded and flux calibrated by mul- 
tiplication with the relative spectral response function 
(RSRF) using the IRS standard star a-Lac for both low 
resolution modules. We built our RSRF by extracting 
the spectra of calibration stars (Cohen et al. 2003) in the 
same way as we perform on our sources, then divide the 
template spectra of those standard stars by the spectra 
we extracted with the column extraction method. Fi- 
nally, we multiplied the spectra of the science target by 
the RSRF, for each nod position. 

Figure 1 shows the positions of the narrower SL and 
wider LL slits for the given observing date overplotted 
on mid-IR images from IRAC, MIPS, and 2MASS in log- 
arithmic scale. Our anticipated slit positions agree quite 
well with the main peak of the mid-IR emission, except 
for NGC 1097, NGC 3310, and NGC 3556, which show 
a more complex morphology. The consequences of slight 
mismatches and extended emission are discussed in the 
following section. 

2.3. Absolute Fluxes and Order Stitching 

After the spectral extraction there was, in some cases, a 
noticeable mismatch between the spectra from the differ- 
ent IRS modules. Since this mismatch is more likely due 
to source flux that was missed in the narrower slits rather 
than unrelated flux that was picked up in the wider slits 
(see below) we scaled the SL2, SL1, and LL2 spectra 
to match the flux density of LL1. The choice of LL1 
as reference slit is appropriate because it has the widest 
slit and largest PSF, and is least sensitive to pointing 
errors or a small spatial extent of the mid-IR emission 
region. We have also used the "bonus orders" SL3 and 
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TABLE 1 
General Properties 



Name 


cx T9nnn a 


^T9nnn a 

J zuuu 


Type 


Refs. b 


D c 


\og(L IR f 


Sl2 C 


S25 c 


Seo e 


Sioo e 








[Mpc] 


[Lq] 


[Jy] 


[Jy] 


[Jy] 


[Jy] 


IC 342 


3 46 48.51 


+68 05 46.0 


SB 


13,24,25 


4.6 


10.17 


14.92 


34.48 


180.80 


391.66 


Mrk 52 


12 25 42.67 


+00 34 20.4 


SB 


2,9,17 


30.1 


10.14 


0.28 


1.05 


4.73 


5.68 


Mrk 266 


13 38 17.69 


+48 16 33.9 


SB+Sy2 


14,20 


115.8 


11.49 


0.32 


1.07 


7.25 


10.11 


NGC 520 


1 24 35.07 


+03 47 32.7 


SB 


4,10,12,13 


30.2 


10.91 


0.90 


3.22 


31.52 


47.37 


NGC 660 


1 43 02.35 


+13 38 44.4 


SB+LINER 


10 


12.3 


10.49 


3.05 


7.30 


65.52 


114.74 


NGC 1097 


2 46 19.08 


-30 16 28.0 


SB+Syl 


19,23 


16.8 


10.71 


2.96 


7.30 


53.35 


104.79 


NGC 1222 


3 08 56.74 


-02 57 18.5 


SB 


1,2 


32.3 


10.60 


0.50 


2.28 


13.06 


15.41 


NGC 1365 


3 33 36.37 


-36 08 25.5 


SB+Sy2 


19,26 


17.9 


11.00 


5.12 


14.28 


94.31 


165.67 


NGC 1614 


4 33 59.85 


-08 34 44.0 


SB 


29,30 


62.6 


11.60 


1.38 


7.50 


32.12 


34.32 


NGC 2146 


6 18 37.71 


+78 21 25.3 


SB 


7,10 


16.5 


11.07 


6.83 


18.81 


146.69 


194.05 


NGC 2623 


8 38 24.08 


+25 45 16.9 


SB 


13,22 


77.4 


11.54 


0.21 


1.81 


23.74 


25.88 


NGC 3256 


10 27 51.27 


-43 54 13.8 


SB 


9,24,25 


35.4 


11.56 


3.57 


15.69 


102.63 


114.31 


NGC 3310 


10 38 45.96 


+53 30 05.3 


SB 


9,10 


19.8 


10.61 


1.54 


5.32 


34.56 


44.19 


NGC 3556 


11 11 30.97 


+55 40 26.8 


SB 


7,10 


13.9 


10.37 


2.29 


4.19 


32.55 


76.90 


NGC 3628 


11 20 17.02 


+13 35 22.2 


SB+LINER 


10,21 


10.0 


10.25 


3.13 


4.85 


54.80 


105.76 


NGC 4088 


12 05 34.19 


+50 32 20.5 


SB 


3,5,10 


13.4 


10.25 


2.06 


3.45 


26.77 


61.68 


NGC 4194 


12 14 09.64 


+54 31 34.6 


SB 


2,9,17 


40.3 


11.06 


0.99 


4.51 


23.20 


25.16 


NGC 4676 


12 46 10.10 


+30 43 55.0 


SB 


15,16 


94.0 


10.88 


0.11 


0.33 


2.67 


5.18 


NGC 4818 


12 56 48.90 


-08 31 31.1 


SB 


5,17 


9.4 


09.75 


0.96 


4.40 


20.12 


26.60 


NGC 4945 


13 05 27.48 


-49 28 05.6 


SB+Sy2 


7,11,31 


3.9 


10.48 


27.74 


42.34 


625.46 


1329.70 


NGC 7252 


22 20 44.77 


-24 40 41.8 


SB 


6,15,18 


66.4 


10.75 


0.24 


0.43 


3.98 


7.02 


NGC 7714 


23 36 14.10 


+02 09 18.6 


SB 


8,27,28 


38.2 


10.72 


0.47 


2.88 


11.16 


12.26 



Note. — Mrk 52 = NGC 4385, Mrk 266 = NGC 5256 
a Commanded coordinates of the slit center, a is given in (h m s), 8 is given in (d m s). 

b References: 1-Ashby et al. (1995), 2-Balzano (1983), 3-Bendo & Joseph (2004), 4-Beswick et al. (2003), 5-Deveraux 
(1989), 6-Fritze-v. Alvensleben & Gerhard (1994), 7-Gao & Solomon (2004), 8-Gonzalcz-Dclgado et al. (1999), 9-Heckman 
et al. (1998), 10-Ho et al. (1997), 11-Iwasawa et al. (1993), 12-Joseph & Wright (1985), 13-Keel (1984), 14-Levenson ct 
al. (2001), 15-Liu & Kennicutt (1995), 16-Lonsdale et al. (1984), 17-Mayya et al. (2004), 18-Miller et al. (1997), 19-Osmer 
et al. (1974), 20-Osterbrock & Dahari (1983), 21-Roberts ct al. (2004), 22-Smith et al. (1998), 23-Storchi-Bergmann ct al. 
(2003), 24-Thornley et al. (2000), 25-Verma et al. (2003), 26-Veron et al. (1980), 27-Weedman ct al. (1981), 28-Brandl et 
al. (2004), 29-Alonso-Herrero et al. (2001), 30-Keto ct al. (1992), 31-Spoon et al. (2000) 

c Distances adopted from Sanders et al. (2003), except for Mrk 52, NGC 4676 and NGC 7252, which were derived from 
measured redshifts via D = 7^(1 + § )> assuming Hq = 71kms~ 1 Mpc -1 . 

d Total 8 — lOOO^im infrared luminosity of the entire galaxy, adopted from Sanders ct al. (2003), except for Mrk 52, NGC 4676 
and NGC 7252, which were derived from measured IRAS fluxes via L IR = 312700 • D 2 ■ 1.8 ■ (13.485i2 Mm + 5.16S25nm + 
2.58560 M m + <Sroo Mm ) where S x is in Jy. 

e IRAS flux densities at 12, 25, 60 and 100/im of the entire galaxy, adopted from Sanders et al. (2003). 



LL3 when they provided better overlap or higher signal- 
to-noise than first and second orders only. The applied 
stitching factors are listed in Table 2 and provide a good 
idea of the uncertainties involved. Since these factors are 
rather large in some cases we would like to emphasize 
the rationale for this approach. The problem of stitch- 
ing together slit apertures of different widths is by no 
means specific to our approach, but applies to basically 
all comparable studies at almost all wavelengths. For 
a black-body-like object with extended, uniform surface 
brightness there will be a jump between the SL slit width 
(3.6") to the LL slit 10.5", corresponding to an increase 
in flux of at least a factor of three. Since all spectral 
energy distributions are continuous, scaling SL to match 
LL seems reasonable to first order. 

However, this approach assumes that the spectral 
properties are not changing within the region covered 
by the LL slit, corresponding to linear scales of about 
200 pc for the nearest objects in our sample the typical 
size of a circum-nuclear starburst over which the spectral 
properties are assumed to not vary substantially. If the 
contributions from a centrally concentrated source, e.g., 
an AGN, were dominant, scaling would lead to an over- 
estimation of the strength of the features originating in 



the center. 

For the latter reason we have decided to list the flux 
densities as measured from the stitched spectra, but not 
to overall scale the IRS spectra to match the spatially 
integrated IRAS flux densities at 25^m. We have calcu- 
lated the vj v average of the IRS spectra over the 25/im 
IRAS filter band. The ratio of IRS to IRAS 25^m flux 
densities is given as the fractional flux 'FF' in Table 2, 
last column. We note that color corrections applied to 
the published IRAS catalog fluxes increase the uncertain- 
ties, but the relative effect on our sample with similar 
SEDs is small. In some cases the ratio is quite small, 
indicating a rather large apparent discrepancy between 
IRS and IRAS. This is mainly due to two reasons: (i) the 
galaxy extends over a large angle and the mid-IR emis- 
sion region is significantly more extended than the IRS 
slit (e.g., NGC 3628, NGC 4945), or (ii) the galaxy has 
significant off-nuclear IR emission peaks (e.g., NGC 3556, 
NGC 4088). 

Figure 1 shows that substructure or extended emis- 
sion on scales of the IRS slit is present in many of our 
objects, most notably in those which have small frac- 
tional fluxes 'FF' in Table 2. The most extreme cases 
are NGC 1097, NGC 3556, and NGC 4088. NGC 3556 
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TABLE 2 
Observational Parameters 



Name 


AOR-kcy 


Obs.date 


*SL a 


tLL* 






Stitching 


factors' 3 






FF C 


IC 342 


9072128 


2004-03-01 


2 X 14 


4x6 


2 


.22 


\ 


.62 


I 


89 





95 





.95 





17 


Mrk 52 


3753216 


2004-01-08 


4x6 


4x6 


1 


.43 


1 


.47 


1 


57 





.95 





.95 





84 


Mrk 266 


3755264 


2004-01-08 


2 x 14 


2 x 14 


1 


.89 


1 


.89 


1 


.79 


o 


.94 


1 


.00 


0. 


58 


NGC 520 


9073408 


2004-07-13 


3 x 14 


3 x 14 


5 


.26 


3 


.53 


4 


11 


o 


95 


1 


.00 





73 


NGC 660 


9070848 


2004-08-07 


4x6 


4x6 


1 


.52 


I 


.33 




45 


o 


.95 


1 


.00 


0. 


73 


NGC 1097 


3758080 


2004-01-08 


2 x 14 


4x6 





.69 





.72 





75 


0. 


69 





.79 


0. 


21 


NGC 1222 


9071872 


2004-08-10 


2 x 14 


2 x 14 


1 


.52 


1 


.18 


1 


.38 


o 


93 





.95 


0. 


76 


NGC 1365 


8767232 


2004-01-04 


4x6 


4x6 


1 


.56 


1 


.09 


1 


34 


o 


88 





.94 





37 


NGC 1614 


3757056 


2004-02-06 


2 x 14 


4x6 


1 


.27 


1 


.02 


1 


19 


o 


88 





.88 


0. 


52 


NGC 2146 


9074432 


2004-02-28 


4x6 


4x6 


2 


.50 


2 


.00 


2 


.28 





.85 





.90 


0. 


55 


NGC 2623 


9072896 


2004-04-19 


2 x 14 


2 x 14 


1 


.25 


1 


.12 


1 


25 


1 


.05 


1 


.05 


0. 


.94 


NGC 3256 


9073920 


2004-05-13 


4x6 


4x6 


1 


.57 


1 


.20 


1 


42 


1 


00 


1 


.00 


0. 


75 


NGC 3310 


9071616 


2004-04-19 


2 x 14 


4x6 





.69 





.76 





76 





.84 





.90 


0. 


28 


NGC 3556 


9070592 


2004-04-18 


2 x 14 


2 x 14 





.85 





.55 





90 





,75 





.82 


0. 


12 


NGC 3628 


9070080 


2004-05-13 


4 x 6 d 


4x6 


1 


.72 


1 


.38 


1 


.66 


0. 


.95 


1 


.00 


0. 


38 


NGC 4088 


9070336 


2004-04-19 


2 x 14 


2 x 14 


1 


.25 


1 


.06 


1 


16 





85 





.90 





13 


NGC 4194 


3757824 


2004-01-08 


2 x 14 


4 x 6 C 


1 


.33 


1 


.20 


1 


.28 





93 





.98 


0. 


81 


NGC 4676 


9073152 


2004-05-15 


6 x 14 


4 x 30 


1 


.30 


1 


.08 


1 


23 





.95 





.97 


0. 


69 


NGC 4818 


9071104 


2004-07-14 


2 x 14 


4 x 6 C 


1 


.05 





.95 


1 


.03 





98 





.98 


0. 


.81 


NGC 4945 


8769280 


2004-03-01 


4x6 


4x6 


1 


.39 





.83 


1 


11 





94 


1 


.03 





28 


NGC 7252 


9074688 


2004-05-15 


3 x 14 


2 x 30 


1 


.54 


1 


.18 


1 


.38 





95 





.98 


0. 


90 


NGC 7714 


3756800 


2003-12-16 


2 x 14 


2 x 14 


1 


.16 





.99 


1 


.10 





.98 


1 


.00 





.78 



a Exposurc time in cycles times seconds. Each cycle in staring mode corresponds to two exposures. Hence, 
to derive the total exposure time one needs to multiply the above numbers by two. 

b Multiplicative factors for the SL1, SL2, SL3, LL2, and LL3 modules, respectively, to stitch the spectral 
fragments together with respect to LL1. The unusually large stitching factors used for NGC 520 are likely due 
to its very irregular and extended structure, which led to a substantial flux loss for the narrowest slits. 

c The fractional flux FF is the ratio of 25/^m flux detected within the LL slit to the total flux of the entire 
galaxy measured by IRAS (Sanders et al. 2003). 

d SL2 was 2 X 14s 

C LL2 was 2 X 14s 



shows a bright, off- nuclear mid-IR source and several IR- 
bright knots along the disk. NGC 4088 is very extended 
with IR emission in the disk that is picked up in the large 
IRAS beam. In the case of NGC 1097 it is clear from Fig- 
ure 1 that the very symmetrical ring is the reason why 
we only see 21% of the flux. Although we do not scale 
our measured spectra to match the IRAS fluxes, the fac- 
tor FF will become very important in section 4 where we 
use absolute fluxes to derive total luminosities and star 
formation rates. 

All of the individual IRS spectra are shown in Fig. 2. 
We have not attempted to correct for the periodic "fring- 
ing" in the spectra longward of about 22/xm, which can 
be very prominent, as in NGC 1614. These artifacts have 
no effect on the analysis carried out in this paper. 

3. ANALYSIS 

Figure 3 shows a normalized overlay of nine starburst 
spectra from our sample with the most prominant spec- 
tral features labeled. The figure illustrates the spectral 
richness of the 5— 38/xm wavelength range and reveals dis- 
tinct differences between individual starbursts. Impor- 
tant common features include PAH emission bands, sili- 
cate emission or absorption features, and emission lines, 
in addition to the information contained in the slope of 
the spectral continuum. In the following subsections we 
discuss how the quantities relevant to our discussion have 
been measured from our spectra. 

3.1. Continuum Fluxes 



In order to characterize the basic properties of the 
spectral continuum we have derived the flux densities 
for three rest frame wavelength ranges: at 5.9 — 6.1/zm, 
at 14.75 — 15.25^m, and at 29.5 — 30.5^m. These wave- 
lengths were chosen to cover a large baseline in wave- 
length while being least affected by PAH emission fea- 
tures, silicate absorption or strong emission lines. The 
flux densities derived for 6/im, 15^im, and 30/zm are the 
median values in the above wavelength ranges, respec- 
tively, and are arguably the best direct estimate of the 
spectral continuum. The measured fluxes are listed in 
Table 3. In section 4.2 we will use these continuum fluxes 
to estimate the total luminosity of the starburst galaxy. 

3.2. Poly cyclic Aromatic Hydrocarbons (PAHs) 

The spectra in Fig. 3 show that the spectral contin- 
uum shape of starburst spectra is dominated by strong 
emission features from PAHs (as previously noted by ISO 
authors, e.g., Genzel & Cesarsky (2000) and references 
therein). Although the first detections of PAHs date 
back in the 1970s, it took more than 10 years to identify 
them. Here we concentrate our analysis on several of the 
strongest PAH features at 6.2/im, 7.7/im, 8.6fim, 11.3/xm, 
14.2/im, and around 17/xm (which is in fact two blended 
PAH complexes centered at 16.4/zm and 17.1/im). Fur- 
ther PAH features can be seen in the spectra (see sec- 
tions 3.5 and 4.6) but have either low S/N or are, at the 
low resolution of the IRS SL+LL modules, blended with 
other features and were thus excluded from this analysis. 
At the IRS low-resolution the 12.7um PAH blends with 
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Fig. 1. — The positions of the SL and LL slits for the given observing date overplotted on IRAC 8^m, MIPS 2Afim or 2MASS K-band 
images for all galaxies in our sample. The scale is logarithmic and overemphasizes the real flux distribution. For discussion see section 2.3. 



the strong [Ne II] line at 12.81/zm. This feature will be 
discussed by Devost et al. (2006) in the context of the 
IRS high resolution spectra. 

The strengths of the 6.2^im, 7.7^m, 8.6/im, and 11.3/zm 
PAH emission bands were derived by integrating the flux 
of the feature in the mean spectra of both nod positions 
above an adopted continuum. For the 6.2 and 11.3/Ltm 
features this baseline was determined by fitting a spline 
function to four or five data points. The wavelength lim- 
its for the integration of the features were approximately 
between 5.94 to 6.56/im in the case of the 6.2/xm PAH 
and between 10.82 and 11.80/xm for the 11.2/um PAH. For 
most spectra this method produced results, reproducable 
to within 5% for repeated fits with different choices of the 
continuum or integration limits to account for uncertain- 
ties within the fitting procedure. Only for the noisiest 
spectra this uncertainty increased to 15%. The baseline 
for the the 7.7/xm and 8.6^m features was derived by fit- 



ting a spline through six data points avoiding small fea- 
tures in the range between 5.5 and 10/im. Our method 
is illustrated in Figure 4. The selected data points for 
the baseline were chosen at the same wavelengths for all 
spectra. Following Peeters et al. (2002) we included one 
point close to 8.2/xm to separate the contribution from 
both features. 

The PAH features at 14.2^m (Hony et al. 2001) and 
17/Ltm (van Kcrckhovcn ct al. 2000) arc intrinsically 
weaker but are located in spectrally less complex re- 
gions, allowing a different approach. The strength of 
the 14.2/im feature was determined through a first or- 
der baseline fit to the continuum at 13.88 — 14.03yum 
and 14.54 — 14.74^m. Similarly, the strength of the 
17/Km complex was determined through a first order 
baseline fit to the continuum at 15.95 — 16.15/im and 
18.20 — 18.40/um. We emphasize that the feature we re- 
fer to as the 17^m PAH is in fact a blend of two PAH 
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Fig. 3.— Overlay of nine IRS starburst spectra of Mkn 52, NGC 520, NGC 1222, NGC 1365, NGC 3310, NGC 3628, NGC 4945, NGC 
7252, and NGC 7714. The objects shown here have been selected to illustrate the full spectral diversity of starbursts with a small number 
of objects. All spectra have been normalized to a flux density of one at 15/^m (dotted vertical line). 



TABLE 3 

Continuum fluxes and extinction 



Name 




FlBum 




T9.8 




[Jy] 


[Jy] 


[Jy] 




IC 342 


0.38 


2.71 


27.26 


0.004 


Mrk 52 


0.03 


0.28 


1.19 


-0.003 


Mrk 266 


0.03 


0.17 


1.16 


0.373 


NGC 520 


0.15 


0.54 


6.12 


0.994 


NGC 660 


0.34 


1.27 


11.52 


1.293 


NGC 1097 


0.12 


0.29 


3.23 


0.130 


NGC 1222 


0.06 


0.38 


3.03 


0.213 


NGC 1365 


0.32 


1.55 


9.61 


-0.034 


NGC 1614 


0.14 


1.10 


5.80 


0.279 


NGC 2146 


0.74 


2.00 


23.07 


0.845 


NGC 2623 


0.06 


0.32 


4.75 


1.544 


NGC 3256 


0.34 


2.36 


21.61 


0.000 


NGC 3310 


0.08 


0.24 


2.79 


0.057 


NGC 3556 


0.02 


0.10 


1.10 


0.230 


NGC 3628 


0.19 


0.46 


4.89 


1.640 


NGC 4088 


0.04 


0.12 


0.76 


0.307 


NGC 4194 


0.22 


0.79 


6.77 


0.371 


NGC 4676 


0.03 


0.07 


0.45 


0.580 


NGC 4818 


0.13 


0.92 


5.56 


-0.160 


NGC 4945 


0.48 


1.83 


49.56 


4.684 


NGC 7252 


0.06 


0.11 


0.69 


-0.098 


NGC 7714 


0.07 


0.53 


3.48 


-0.096 



complexes at 16.4l«ti and 17.lL!m. The latter is furthc- 
more contaminated by the H2S(1) line at 17.0/im. The 
individual components of this complex cannot be prop- 
erly resolved in the low-resolution spectra and we give 
only combined fluxes here. The equivalent widths for all 
features were derived by dividing the integrated PAH flux 
above the adopted continuum by the continuum flux den- 
sity at the center of the feature (indicated by the vertical 
dot-dashed line in Figure 4) . The derived PAH fluxes and 
equivalent widths arc listed in Table 4. 

It is important to note that the values in Table 4 have 
not been corrected for extinction. While in many sources 
the silicate absorption band around 9.8/im is very weak, 
objects severely affected by extinction, like NGC 4945, 
show a strong 7.7^m PAH but much weaker 8.6/xm and 
11.3y«m features. The cause of these variations will be 
addressed in section 4.6. 

We do not list uncertainties for the PAH strengths in 
Table 4. The statistical errors from the fits to the high 
S/N are small compared to other uncertainties such as: 
(i) The error in the absolute flux calibration of each mod- 
ule, which is currently about 5 — 10%. (ii) The error from 
scaling the individual orders for an extended source to 
match, as discussed in section 2.3. (Hi) The error in 
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TABLE 4 

Main PAH feature strengths shortward of 18/itM 



6.2/un PAH 7.7/im PAH 8.6/im PAH 11.3/mi PAH 14.2/im PAH 17//m PAH 
Name Flux a EW b Flux a EW b Flux a EW b Flux a EW b Flux a EW b Flux a EW b 



IC342 


14.65 


0.497 


31.75 


0.581 


7.39 


0.168 


18.36 


0.492 


0.50 


0.014 


9.70 


0.178 


Mrk52 


1.07 


0.535 


2.14 


0.552 


0.51 


0.151 


1.02 


0.316 


0.02 


0.005 


0.66 


0.142 


MrkioD 


0.59 


0.619 


0.92 


0.467 


0.00 


0.126 


0.47 


0.422 


0.08 


0.065 


0.51 


0.403 


NGC0520 


5.60 


0.563 


13.97 


0.528 


1.99 


0.126 


5.62 


0.798 


0.20 


0.024 


2.91 


0.494 


NGC0660 


12.91 


0.504 


27.36 


0.518 


3.81 


0.123 


9.72 


0.701 


0.43 


0.023 


6.26 


0.365 


NGC1097 


5.04 


0.459 


9.30 


0.488 


2.26 


0.168 


5.34 


0.661 


0.12 


0.021 


3.82 


0.523 


Mnm nan 

NGC1222 


2.17 


0.624 


4.64 


0.606 


0.79 


0.130 


2.56 


0.566 


0.08 


0.015 


1.47 


0.231 


NGC1365 


2.74 


0.111 


6.37 


0.213 


1.30 


0.045 


4.60 


0.180 


0.62 


0.029 


5.59 


0.250 


NGC1614 


13.05 


0.561 


24.67 


0.514 


5.01 


0.138 


9.77 


0.379 






5.10 


0.154 


NGC2146 


22.22 


0.545 


58.26 


0.643 


10.51 


0.175 


21.53 


0.829 


1.20 


0.044 


13.60 


0.492 


NGC2623 


1.83 


0.598 


4.00 


0.454 


0.76 


0.131 


1.28 


0.527 


0.01 


0.002 


0.71 


0.205 


NGC3256 


8.02 


0.603 


17.60 


0.533 


3.03 


0.123 


8.52 


0.471 


0.71 


0.040 


8.12 


0.205 


NGC3310 


3.35 


0.789 


5.38 


0.591 


1.16 


0.178 


2.98 


0.748 


0.07 


0.022 


1.08 


0.229 


NGC3556 


1.22 


0.502 


2.92 


0.523 


0.51 


0.135 


1.52 


0.811 


0.06 


0.039 


0.84 


0.542 


NGC3628 


7.45 


0.500 


19.72 


0.588 


1.59 


0.095 


4.14 


0.797 


0.25 


0.034 


3.45 


0.684 


NGC4088 


1.17 


0.496 


2.43 


0.483 


0.51 


0.130 


1.27 


0.603 


0.05 


0.025 


0.97 


0.509 


NGC4194 


7.09 


0.529 


14.67 


0.578 


3.06 


0.165 


6.31 


0.590 


0.25 


0.022 


3.10 


0.241 


NGC4676 


1.05 


0.610 


2.15 


0.551 


0.44 


0.192 


0.97 


0.812 


0.04 


0.038 


0.52 


0.590 


NGC4818 


4.06 


0.459 


9.56 


0.555 


1.84 


0.123 


4.39 


0.344 


0.05 


0.004 


2.43 


0.151 


NGC4945 


12.13 


0.432 


40.34 


0.490 


0.13 


0.003 


3.55 


0.558 


0.50 


0.024 


5.37 


0.519 


NGC7252 


2.09 


0.585 


4.48 


0.549 


1.06 


0.176 


2.95 


0.931 


0.04 


0.024 


1.35 


0.792 


NGC7714 


2.67 


0.601 


5.62 


0.642 


1.06 


0.135 


2.79 


0.394 


0.04 


0.006 


1.05 


0.114 



a Flux in units of 10 19 Wcm 2 
b Equivalent width in units of /im 



3.3. Emission Lines 

The 5 — 38^m wavelength range contains numerous 
strong emission lines. Among those are the following 
forbidden lines - sorted by wavelength and with their 
excitation potentials in parentheses: 
[Ar II] 6.99^m, [At III] 8.99^m, [S IV] 10.51/rai, 
[Nell] 12.81/xm, [Nelll] 15.56/zm, [S III] 18.71//m, 
[OTV] 25.89/xm, [S III] 33.48/zm, and [Sill] 34.82/zm. 
In addition, we detect the pure rotational lines of 
molecular hydrogen H2 (0,0) S(5) 6.91/im (blended with 
[Aril]), H 2 (0,0) S(3) 9.66^m, H 2 (0,0) S(2) 12.28/xm, and 
H 2 (0,0) S(l) 17.03^m. 

All of these lines have been identified and labeled in 
Fig. 3. We list them here since they can be easily de- 
tected, even at R < 100. However, the flux measure- 
ments of the fine-structure lines can be done much more 
accurately from the IRS high-resolution spectra, which 
is the subject of a complementary paper discussing the 
ionic properties of the ISM (Devost et al. 2006). 

3.4. Silicate Absorption and Optical Depth 

The wavelength coverage of IRS is ideally suited for 
a detailed study of the strong vibrational resonances in 
the silicate mineral component of interstellar dust grains. 
Amorphous silicates — the most common form of sili- 
cates — have a broad Si-0 stretching resonance peaking 
at 9.8/um and an even broader O-Si-0 bending mode 
resonance peaking at 18.5/im. 

We have estimated the apparent optical depth in the 
9.8/im silicate feature from the ratio of the local mid- 
infrared continuum to the observed flux at 9.8/im. For 
shallow silicate features, the local continuum may be de- 
fined as an Fa power law interpolation between contin- 
uum pivots at 5.5/xm (averaged 5.3 — 5.7/im flux) and 
14.5/um (averaged 14.0 — 15.0/im flux), thus avoiding the 




Rest Wavelength [fj.m] 



Fig. 4. — Illustration of how the 7.7/im PAH strength was mea- 
sured in NGC 7714. The light grey dots represent the continuum 
points, interpolated by the dashed spline function. The vertical 
dotted lines indicate the integration range for the 7.7/^m feature. 



the approximation of the underlying continuum, which 
is often dominated by strong, adjacent emission and ab- 
sorption features and varies from source to source, and 
(iv) The error in the strength of the 17/im PAH from the 
blending of two PAH complexes and the H 2 S(1) line. 

Our measurement procedures have been designed to 
minimize these errors as best as possible. The dominant 
error remains the uncertainty in the flux calibration, and 
we estimate the error in the PAH measurements to be in 
the order of 10%. 
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TABLE 5 
Relative absorption values 
interpolated from draine (1989) and 
Lutz (1999) 



Wavelength 


Af raine /A v 


A\ utz IAy 


6.2/^m 


0.0164 


0.0489 


7.7jum 


0.0112 


0.0440 


9.8/jm 


0.0554 


0.1289 


11.3//m 


0.0375 


0.0872 



main PAH emission complexes at 6 — 9/iin and 11 — 13yum. 

For starburst spectra with a more pronounced silicate 
feature, the continuum in the 14— 15/im range is affected 
by weak absorption from the overlapping wings of the 9.8 
and 18/im silicate features. For these spectra we define 
a second local continuum, by replacing the continuum 
pivot at 14.5/zm by a continuum pivot at 28/xm (averaged 
27.5 — 28.5^m flux) and use the average of the optical 
depths derived with either local continuum as the best 
estimate of the apparent 9.8/im silicate optical depth. 

Most fluxes used in this analysis are observed flux den- 
sities. However, in Fig. 15 we correct the observed fluxes 
for extinction. To get an estimate of the uncertainties 
we use two extinction laws from Draine (1989) and Lutz 
(1999) and show the difference in Fig. 15. The relative 
silicate absorption values A\jAy for the relevant PAH 
wavelengths are listed in Table 5. 

Some estimates of Tg. 8Mm in Table 3 are negative, im- 
plying that silicates are observed in emission. However, 
the absolute values are quite small and may just rep- 
resent uncertainties in our baseline definition (see sec- 
tion 4.3 for a discussion). We also note that the true 
silicate optical depth may be significantly larger than 
the apparent silicate optical depth if the emitting and 
absorbing sources are mixed along the line of sight; if 
part of the silicate column is warm; or if the absorption 
spectrum is diluted by unrelated foreground emission. 

3.5. Spectral Features in the 5 — 8/xm Range 

The 5 — 8/xm spectral range of starburst galaxies is 
extremely rich in atomic and molecular emission and ab- 
sorption features, and dominated by emission from the 
6.2/im PAH feature and the blue wing of the 7.7/im PAH 
complex. Weaker emission features are expected from 
atomic lines ([Fell] at 5.34/tm and [Aril] at 6.99/im), 
molecular hydrogen (H 2 S(7) at 5.51/xm and H 2 S(5) at 
6.91/xm), and 'combination-mode' PAH emission bands 
(at 5.25/im and 5.70/im). Absorption features of wa- 
ter ice and hydrocarbons, commonly detected in deeply 
obscured galactic nuclei (Spoon et al. 2002), would be 
expected at 6.0/^m (water ice) and 6.90/mi and 7.25/im 
(C-H bending modes in aliphatic hydrocarbons). 

As illustrated by the average starburst spectrum in 
Fig. 6, PAH combination- mode emission features are 
common in our starburst spectra. Inspection of the indi- 
vidual spectra shows that the 5.25/xm feature is usually 
double-peaked due to blending with the 5.34/xm [Fell] 
line (e.g. NGC 1222 and NGC 3256). Likewise, the 
profile of the 5.70/im PAH feature is affected by the 
presence of the H 2 S(7) line at 5.51/um (most notably 



NGC 2623 and NGC 4945). The 5.70/xm PAH feature 
appears strongest in the spectrum of NGC 4945 (Fig. 5). 
The ratio of the 5.7[im to 6.2/im PAH in this source is 
0.29, about 5 times higher than for most other starburst 
galaxies in our sample. Interestingly, the red wing of the 
5.7 fim PAH feature coincides with the steep onset of the 
6/zm water ice absorption feature, as illustrated in Fig. 5 
by the steep change in slope at 5.7/im in the spectrum 
of the ULIRG IRAS 20100-4156. Simple spectral mod- 
eling confirms that a screen of water ice absorption can 
indeed mimick a stronger 5.7/nm PAH feature by sup- 
pressing the red wing of the feature and the adjacent 
5.9/im continuum. The presence of water ice in the nu- 
cleus of NGC 4945 is further supported by the discovery 
of a 3fim water ice absorption feature in the ISO-PHT 
S and VLT-ISAAC spectra of the nucleus (Spoon et al. 
2000). For the remaining galaxies in our sample, the 
5 — 6/um spectral structure does not provide evidence for 
the presence of water ice absorption. Hence, apart from 
NGC 4945, shielded cold molecular gas may not be as 
abundant in starburst nuclei as in ULIRG nuclei. 
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Fig. 5. — Comparison of the average starburst template (dotted 
line) from Fig. 6 to the most extincted source within our sam- 
ple, NGC 4945 (dashed line), and to the ULIRG IRAS 20100-4156 
(solid line). The spectra are normalized at 5.7/im. IRAS 20100- 
4156 has been chosen as a representative ULIRG from the sample 
of Spoon et al. (2006). 

Absorption features of aliphatic hydrocarbons at 6.85 
and 7.25/xm are thought to be tracers of the diffuse ISM 
(Chiar et al. 2000). These features are readily detected 
in the spectra of deeply obscured ULIRG nuclei such as 
IRAS 20100-4156 (Fig. 5). The average starburst spec- 
trum (Fig. 6), in contrast, does not show similarly pro- 
nounced structure and the 7/zm range is dominated in- 
stead by the blend of 6.91/um H 2 S(5) and 6.99/um [Aril]. 
However, individual starburst spectra show weak spec- 
tral structure in the 6.5 — 7.0yum range. At 6.65/um, 
a weak emission feature seems to be present, most no- 
tably in the spectra of IC 342, NGC 660, NGC 1614, 
NGC 2146, NGC 4088, NGC 4945 and NGC 7252. 
Pecters et al. (2002) identify this feature as a PAH emis- 
sion band. The red wing of the 6.65/zm emission feature 
lies close to the expected onset of the 6.85/xm hydrocar- 
bon absorption feature (Fig. 5). Especially in the spec- 
trum of NGC 4945, the 6.5 — 7.0/xm spectral structure 
is consistent with the presence of hydrocarbon absorp- 
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tion at a strength of r(6.85/xm)=0.15±0.05 (Fig. 5). For 
other starburst galaxies the spectral structure is too shal- 
low and/or the S/N of the spectra too low to identify 
hydrocarbon absorption with sufficient confidence. 



[Aril] [Arlll] [S IV] [Ne II] [Ne III] [S III] [O IV] 




5 10 20 30 
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Fig. 6. — Average IRS spectrum of 13 starburst galaxies (IC 
342, NGC 660, NGC 1097, NGC 1222, NGC 2146, NGC 3310, 
NGC 3556, NGC 4088, NGC 4194, NGC 4676, NGC 4818, NGC 
7252, and NGC 7714). All spectra have been normalized to a flux 
density of one at 15/im before co-addition. The dotted line shows 
the ISO-SWS spectrum of M82 (Sturm et al. 2000) for comparison. 



3.6. The Starburst "Template Spectrum" 

Although many spectral features show variations from 
one starburst galaxy to another, the magnitude of 
these changes is relatively small compared to the differ- 
ences between different classes of objects, such as AGN, 
quasars, ULIRGs, or normal galaxies. It is common prac- 
tice to classify objects in these categories, and the avail- 
ability of reference spectra is of great interest. A refer- 
ence "template" spectrum would allow a comparison how 
close a given spectrum is to a typical starburst galaxy or 
if it shows any atypical features. Furthermore, classifi- 
cations of objects at high redshift often require fits to a 
library of template spectra. 

We have constructed a high signal-to-noise "starburst 
template" from the spectra of IC 342, NGC 660, NGC 
1097, NGC 1222, NGC 2146, NGC 3310, NGC 3556, 
NGC 4088, NGC 4194, NGC 4676, NGC 4818, NGC 
7252, and NGC 7714. These are basically the objects 
from our sample with high fluxes and without a strong 
AGN component. The spectra have been normalized to a 
flux density of unity at 15/im before averaging. Figure 6 
shows the resulting "template spectrum" and Table 6 
lists the "average" spectral properties derived from this 
composed spectrum. 

Figure 6 also shows the ISO-SWS spectrum of M82 
(Sturm et al. 2000), which is often being used as a star- 
burst template, for comparison. Although the spectral 
slope longward of 15/im is very similar, the two spectra 
show several distinct differences: the ISO-SWS spectrum 
of M82 docs not show the pronounced PAH complex 
around 17/im, it shows much stronger silicate absorp- 
tion, and the flux density shortward of 12/im is almost 
a factor of two lower that in our average starburst tem- 
plate. We provide the spectrum in ascii table format on 



TABLE 6 
Properties of the "average" 
starburst galaxy 



Parameter 


Value 




26% 


Fl5/J.m 


100% 




856% 


6.2^m PAH EW 


0.53/im 


7.7^m PAH EW 


0.53/im 


8.6^m PAH EW 


0.15/im 


11.3/wn PAH EW 


0.66/im 


14.2/mi PAH EW 


0.02^m 


17/im PAH EW 


0.36/im 


T9.8 


0.24 ±0.10 



our website 5 . 

4. RESULTS AND DISCUSSION 

It is important to keep in mind that the starburst 
spectra presented in this paper represent an entire star- 
burst region, including numerous (super-)star clusters at 
various ages and evolutionary states, the surrounding 
PDRs which are internally and externally excited, and 
the warm and cold dust spread across the entire region as 
well as localized dust condensations. While many of the 
properties of local substructures are expected to vary sig- 
nificantly, the overall significance of these variations may 
be averaged out in the spatially integrated spectra. Our 
aim here is to search for global trends between the spec- 
tral properties derived in the previous section (silicate 
absorption features, PAH features, spectral continuum) 
and the global starburst properties (Ljr, radiation field). 

4.1. The Continuum Slope as a SB/ AGN Diagnostic 

The slope of the mid-IR spectral continuum depends 
on the optical thickness, composition and temperature 
of graphite dust grains, which are related to the amount 
of silicate grains (Mathis, Rumpl & Nordsieck 1977). 
The dominating species in the IRS spectral range are hot 
(> 100 K), large grains heated by ionizing, non-ionizing, 
and Lya photons inside Hil regions, and small (< 100 A) 
grains heated by non-ionizing photons outside Hn re- 
gion (Mouri et al. 1997). In Fig. 7 we plot the contin- 
uum slope, parametrized by the ratio of the 15/^m and 
30/im flux densities versus the optical depth at 9.8/im 
(section 3.4). The filled symbols in Fig. 7 (and all other 
figures thereafter) correspond to starbursts with a weak 
AGN component as identified in Table 1. Although not a 
tight correlation one can see the general trend that star- 
bursts with stronger silicate absorption tend to have a 
steeper continuum at longer wavelengths. 

The slope of the dust continuum depends on the en- 
ergy distribution and spatial concentration of the heat- 
ing source(s). Dale et al. (2000) found from ISOCAM 
data of 61 galaxies a 6.75/im / 15/im continuum slope 
near unity for more quiescent galaxies, whereas that ra- 
tio drops (i.e., the slope steepens) for galaxies with in- 
creased starburst activity. Furthermore, the continuum 
slope can serve as a discriminator between massive stars 

5 http://www.strw.leidenunlv.nl/ brandl/SBtemplate.html 
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Fig. 7. — Optical depth at 9.8/im versus the slope of the con- 
tinuum as measured by the flux ratio of f 5fim to 30^tm. The filled 
diamonds correspond to star bursts with a weak AGN component. 
The dotted line is a linear fit to the data points. The highly ob- 
scured source NGC 4945 is not included in the fit.. 



or an AGN as the underlying power source. This has al- 
ready been known since IRAS (e.g., Wang (1992)), and 
was further refined in numerous papers based on ISO- 
CAM observations. (See Genzel & Cesarsky (2000) for 
a more comprehensive overview). For instance, Laurent 
et al. (2000) studied the ISOCAM colors of a large vari- 
ety of cxtragalactic objects revealing clear general trends 
between AGN, photo-dissociation region (PDR), and Hn 
region dominated spectra. However, the characterization 
of an individual object often turns out to be difficult due 
to the intrinsically large scatter: the IRAS 12/zm filter in- 
cludes the silicate absorption band as well as several PAH 
emission features and strong emission lines. The wave- 
length range covered by ISOCAM is limited to shorter 
wavelengths which are dominated by hot dust and a large 
variety of emission and absorption features (compare to 
Fig. 5). 

Our narrowband continuum fluxes largely avoid these 
problems. Figure 8, analogous to a color-color diagram 
in stellar astronomy, shows the continuum flux ratios at 
6/im/15^m versus 15/zm/30^m. The total infrared lumi- 
nosity Lir is represented by the size of the symbols. For 
comparison the figure also contains the classical AGN 
CenA (Sy2), IZwl (Syl), Mrk 3 (Sy2), NGC 1275 
(Sy2), NGC 4151 (Syl.5), and NGC 7469 (Syl.2) from 
Weedman et al. (2005). We also show the starburst- 
dominated ULIRGs IRAS 12112+0305, IRAS 22491- 
1808, and the AGN-dominated Mrk 231, Mrk 463, and 
Mrk 1014 from Armus et al. (2006b). NGC 6240 is a pe- 
culiar case with an intrinsic fractional AGN contribution 
to the bolometric luminosity of 20 — 24% (Armus et al. 
2006a). 

Several conclusions can be drawn from Fig. 8. First, 
there is a large scatter along the j/-axis with no obvi- 
ous correlation with total starburst luminosity or galaxy 
type. Hence, the 6/zm/15/im continuum flux ratio does 
not appear to be a good diagnostic. Second, using 
the 15/im/30^im continuum flux ratio, classical (strong) 
AGN can be clearly separated from the starburst galaxies 
(including the ones with weak nuclear activity), with the 
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Fig. 8. — "Color-color plot" of the 6^tm/f5/im versus 
15/im/30/im flux densities. The size of the diamonds, representing 
our starburst sample, increases linearly with log L jr. The filled 
diamonds correspond to starbursts with a weak AGN component 
and the square refers to the average starburst template (Table 6). 
For comparison the figure also contains AGN (pluses) from Weed- 
man et al. (2005), as well as ULIRGs (filled circles) from Armus et 
al. (2006a) and Armus et al. (2006b). 



AGN having a significantly shallower mid-IR spectrum. 
The grey-shaded bar in Fig. 8 indicates the transition 
region between AGN and starburst-dominated systems, 
and lies approximately at 

0.17 < i 7 i5 Aim /i 7 30^m < 0.21. 

Third, the difference in spectral slope also seems to apply 
to ULIRGs depending on their dominant power source. 
Hence, this technique may have a much broader applica- 
tion and should be verified with a much larger sample of 
different classes of objects. 

4.2. The Continuum Fluxes as Measures of Lj R 

The total infrared luminosity is an important parame- 
ter to estimate the energetics of a starburst and to char- 
acterize the underlying stellar population and rate of star 
formation. Lm is often derived from the four IRAS filter 
bands (Sanders & Mirabel 1996). In this subsection we 
will check how accurately Lm can be derived from the 
two IRS continuum fluxes at Fi 5Mm and F 30t j,m alone. 

Numerous attempts to extrapolate Ljr from one or 
two, mainly broad band fluxes can be found in the lit- 
erature. Takeuchi et al. (2005) discuss various esti- 
mators of infrared luminosities, and find - for a very 
large sample of 1420 galaxies of different type - correla- 



tions of the form \ogL IR = 1.02 + 0.972 log L 



and 



logi/_R = 2.01 + 0.878 logL 25t i m - Both L 12 ^ m and L 25f im 
allow to predict Lj R to an accuracy within a factor of 
4 — 5 at the 95% confidence level over a wide range in 
luminosities. These uncertainties are similar to the lu- 
minosities derived from the MIPS 24/im flux alone for a 
large sample of SINGS galaxies (Dale et al. 2005). Chary 
& Elbaz (2001) found a similar relation fitting 15^m 
ISO fluxes: \ogL IR = (1.05 ± 0.174) + 0.998 log L 15fim . 
Forster Schreiber et al. (2004) found that the monochro- 
matic 15/im continuum emission is directly proportional 
to the ionizing photon luminosity, and hence the total 
infrared luminosity. 
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As discussed in section 4.1 the IRS fluxes fi5^ m and 
-p30/im provide a rather accurate estimate of the "true" 
spectral continuum. In Figure 9 we plot a combination of 
FiSfim and F 30f mi times D 2 versus the total infrared lu- 
minosity derived from the IRAS bands. Since the IRAS 
beam usually covers the entire starburst while the nar- 
rower IRS slits can only collect a fraction of the total 
luminosity for local starbursts (Fig. 1) we have corrected 
the observed F 15fim and F 30Mm (indicated by the pluses 
in Fig. 9) for the slit losses (see section 2.3 for more de- 
tails). A least-square fit to the corrected fluxes yields: 

Fir = F> 2 ■ (4.27 • Fi 5fml + 11 • F 30Aim ) 
where D is the distance in kiloparsecs, i*i5^ TO and F^^ m 
the IRS flux densities in Janskies. 
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Fig. 9. — The estimated infrared luminosity L||| based on the 
aperture loss corrected IRS continuum fluxes Fi5 Mm and F3o Mm for 
a source at distance D [kpc] versus Lju from IRAS. The dashed 
line represents a linear fit to the data (in log- log space), and the 
fitted function is given in the upper left, together with the stan- 
dard error a. The uncorrected fluxes are indicated by gray plus 
signs. The filled symbols correspond to starbursts with a weak 
AGN component. 



The correlation in Fig. 9 is extremely tight, including 
the weak AGN, with a standard error (mean scatter) of 
only 0.09 in log-space, i.e., the IKS estimated infrared lu- 
minosities agree within 23% with Lir. This is much more 
accurate than the estimates by Chary & Elbaz (2001) 
and Takeuchi et al. (2005). The excellent correlation sug- 
gests that, at least for a homogeneous sample of starburst 
galaxies, Fi 5)lm and F 30flm can be used to accurately de- 
rive L IR . 

4.3. A Large Variety in Silicate Absorption 

Figure 10 shows, from top to bottom, a series of star- 
burst spectra, sorted by increasing absorption of the 
9.8^m silicate resonance. Since the peak of the reso- 
nance coincides with a minimum between the 7 — 9^m 
and 11 — 13/im PAH emission complexes, the effect of 
silicate absorption only becomes apparent towards the 
lower half of the plot. Nevertheless, the figure strikingly 
illustrates the strong effect of amorphous silicates on the 
overall 5 — 38/^m spectral shape and the large variations 
even within one class of objects. Fig. 10 shows that the 
10/xm trough can in fact become a dominating feature of 
the spectral shape of starbursts. 
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Fig. 10.— The IRS spectra of Mkn 52, NGC 7714, NGC 1222, 
NGC 7252, NGC 2146, NGC 520, NGC 3628, and NGC 4945 ar- 
ranged to illustrate the gradual effect of increasing silicate absorp- 
tion from top to bottom. The flux densities F\ have been arbitrar- 
ily scaled for better comparision. 



The top two spectra also show a broad emission struc- 
ture beginning at about 16/xm and extending to 21 — 
25/xm. The shape of this feature is consistent with that 
of an 18/im silicate emission feature, which would make 
these the first detections of this feature in starburst 
galaxies. However, the emission structure could also 
be due to the C-C-C in-plane and out-of-plane bend- 
ing modes of PAHs (van Kcrckhovcn et al. 2000). The 
identification of this feature will be discussed in a future 
paper. The remaining starburst galaxies in Fig. 10 do 
not show evidence for a similar emission structure. The 
13 - 35^m spectra, from NGC 1222 (top) to NGC 4945 
at the bottom, show an increasingly pronounced depres- 
sion, peaking at 18.5/im, signalling increasingly strong 
silicate absorption. The latter result is in full agreement 
with the trend found for the 9.8^m silicate feature. 

Among all the galaxies classified as starbursts, 
NGC 4945 (bottom spectrum in Fig. 10) is a "special 
case" exhibiting by far the strongest dust obscuration to 
its nuclear region (Spoon et al. 2000). Based on the IRS 
spectrum, the apparent optical depth in the 9.8^m sili- 
cate feature is at least four and may be higher, depending 
on the choice of the local continuum. Apart from strong 
amorphous silicate absorption, the line of sight also re- 
veals the presence of a 23^tm absorption feature. Follow- 
ing the analysis of deeply obscured lines of sight towards 
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ULIRG nuclei (Spoon et al. 2006), we attribute the 23/tm 
feature to crystalline silicates (forsteritc). Their detec- 
tion in NGC 4945 suggests that crystalline silicates are 
perhaps a more common component of the ISM and not 
just limited to ULIRG nuclei. 

Figure 11 compares the average starburst template 
from Fig. 6 to NGC 4945, the most extincted source 
within our sample. For comparison we also show the 
heavily embedded ULIRG IRAS 08572+3915 from Spoon 
et al. (2006). The usually rather shallow 18/xm sil- 
icate band reduces the continuum by a large factor 
compared to the average starburst spectrum. While 
NGC 4945 has a similarly strong silicate absorption fea- 
ture as the ULIRG IRAS 08572+3915, the starburst 
also shows strong PAH complexes near 16.4/tm and 
17.1/tm, which are absent in the ULIRG spectrum. Crys- 
talline silicates, causing the features at 16 and 19/im in 
IRAS 08572+3915, appear to be absent in the starburst 
template, although they are more difficult to discern at 
our spectral resolution given the presence of PAH fea- 
tures and the [S III] line in that range. 



region. Table 4 already indicates that the relative fluxes 
of individual PAH features are not constant for different 
starbursts. Lu et al. (2003) have found a 25% spread in 
the ratio of the 11.3/im/ 7.7/mi PAH fluxes which they 
attribute to intrinsic galaxy-to-galaxy variations. How- 
ever, extinction can affect the relative strength of fea- 
tures at different wavelengths by different amounts (Ta- 
ble 5). Rigopoulou et al. (1999) found that the domi- 
nant influence on the (ULIRG) PAH ratio is extinction, 
and that the 6.2/im PAH gets suppressed relative to the 
7.7/tm PAH for heavily dust extincted systems. In par- 
ticular, the 8.6/xm and 11.3/im PAH features lie at wave- 
lengths that are heavily affected by silicate absorption. 

Fig. 12 shows the 11.3/im / 7.7/im and 11.3/im / 6.2/im 
PAH flux ratios versus the apparent silicate optical depth 
Tg.8- A trend that starbursts with stronger dust ab- 
sorption show relatively weaker 11.3/im PAHs is evident. 
Fig. 12 suggests that extinction can affect the relative 
PAH strength in starbursts by up to a factor of about 
two. 
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Fig. 11. — Comparison of the average starburst template (dotted 
line) from Fig. 6 to the most extincted source within our sample, 
NGC 4945 (dashed line), and to the ULIRG IRAS 08572+3915 
(solid line), chosen as an extreme case of extinction from the sam- 
ple of Spoon et al. (2006). The spectra have been normalized at 
14.5/iin. 

While NGC 4945 shows by far the most extreme sil- 
icate absorption within our sample it is not an exotic 
object but defines the endpoint of a sequence of increas- 
ing optical depths. Recently, Dale et al. (2006) reported 
on the Spitzer-SINGS survey of 75 nearby galaxies. The 
typical target in their sample shows only modest dust ob- 
scuration, consistent with an Ay ~ 1 foreground screen 
and a lack of dense clumps of highly obscured gas and 
dust. Unfortunately, we cannot distinguish between the 
local and global dust distributions. The differences may 
be due to evolutionary states, local geometries or other 
reasons. However, it is evident that starburst galaxies 
can have very little or large amounts of extinction, and 
the presence/absence of a strong dust feature is not a 
characteristic item. 

4.4. The Mixture of PAHs and Dust 

In this section we will investigate how the derived PAH 
strengths may depend on extinction within the starburst 
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Fig. 12. — The PAH flux ratios 11.3/im/7.7/um (diamonds) and 
11.3/tm/6.2jum (filled circles, divided by a factor 2.05 to match the 
other ratio) versus the optical depth at 9.8/im, rg.g from Table 3. 
The dashed line shows a linear fit to both ratios but excluding 
NGC 4945. 



In the remainder of this subsection we will investigate if 
the PAH equivalent width is related to the total infrared 
luminosity of the starburst. Rigopoulou et al. (1999) 
found a ratio of L pah / Lm which is similar for starburst- 
dominatcd ULIRGs and for template starbursts, i.e., no 
dependency on Ljp. In contrast, Lu et al. (2003) found 
a steady decrease of the PAH strength with increasing 
IR- activity. Figure 13 shows the equivalent widths of 
both the 6.2/im and 7.7/im PAH as a function of the to- 
tal infrared luminosity for our starburst sample. Within 
the uncertainties, indicated by the scatter between the 
6.2/xm and 7.7/tm PAHs for the same object, the PAH 
equivalent widths remain constant over a factor of 50 in 
total luminosity. This finding is in good agreement with 
Peeters et al. (2004) (and references therein) who found 
that the fraction of the total PAH flux emitted in the 
6.2/im PAH band varies only slightly with an average of 
28% ± 4%. In other words, the PAH flux and the un- 
derlying warm dust continuum scale proportionally, and 
the PAHs and dust must be well mixed, at least on large 
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scales, to show these correlations. This rules out a sce- 
nario of a luminous, dusty nucleus surrounded by a large 
PDR with little extinction, in favor of smaller, clumpier 
structures. 
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Fig. 13. — The 6.2/j.m (diamonds) and 7.7/im (filled circles) PAH 
equivalent widths versus infrared luminosity Ljj^. The dashed, 
horizontal line indicates a zeroth order polynomial fit to all data 
points, excluding the outlier NGC 1365. The filled diamonds cor- 
respond to starbursts with an AGN component. 



4.5. PAH Luminosity and Star Formation Rate 

Kennicutt (1998) has shown that the 8 — 1000/xm in- 
frared luminosity Lir of starbursts is a good measure of 
the star formation rate (SFR) given by 

SFR[M Q yr _1 ] = 4.5 x KT^L/^erg s" 1 ]. 

The SFR determines the number of young, massive stars, 
which provide the (far-)UV photons to excite both PAH 
molecules and dust grains. If both species get excited 
by the same photons PAHs could potentially be used as 
quantitative tracers of the SFR (see section 4.2 for the 
dust luminosity). 

Generally, PAHs arc considered the most efficient 
species for photoelectric heating (Bakes & Tielens 1994) 
in the PDRs. Molecular gas in these boundary layers, 
surrounding the Hn regions, is exposed to far-UV radi- 
ation (6 — 13.6 eV), which strongly influences its chemi- 
cal and thermal structure (Tielens & Hollenbach 1985). 
PAHs are stochastically heated by these UV photons, 
predominantly originating from massive stars, and hence 
expected to be good tracers of star formation. On large 
angular scales - similar to our case - Forster Schreiber 
et al. (2004) found that the 5 — 8.5/xm PAH emission 
constitutes an excellent indicator of the star formation 
rate in circumnuclear regions and starbursts as quanti- 
fied by the Lyman continum flux, i.e. in regions where 
the energy output is dominated by massive star forma- 
tion. However, it has been known for a long time that 
PAHS can also be excited by visible photons (e.g., Uchida 
et al. (1998)), and that PAHs can trace other sources 
besides massive young stars, such as planetary nebulae 
and reflection nebulae. If the observed PAH flux is in- 
tegrated over the whole galaxy it may predominantly 
trace B stars, which dominate the Galactic stellar energy 



budget, rather than very recent massive star formation 
(Peeters et al. 2004). This is in agreement with higher 
angular resolution observations of the 3.3^m PAH fea- 
ture at the VLT by Tacconi-Garman et al. (2005), who 
found a decrease in the PAH/continuum ratio at the sites 
of the most recent star formation. 
In section 4.4 we have seen that the PAH equivalent 
~ z width does not dependent on Lir. However, that does 
: not necessarily mean that PAHs aren't good quantitative 
j tracers of star formation if the PAH feature and the un- 
: derlying continuum scale proportionally. In Fig. 14 we 
compare the flux in the 6.2/im PAH feature against the 
total infrared luminosity Lir from IRAS. To search for 
I a physically meaningful correlation one needs to take the 
: distance of the object into account as well as the fact that 
the narrow IRS slit misses some of the total flux. Hence 
: we multiply the measured PAH fluxes with the square of 
J the distances and divide by the fractional flux factor FF 
(Table 2). A remarkably good fit can been achieved with 

j log (Lfi H ) = 1.13 x log (F 6 . 2 ^ mPAH D 2 ) 

120 where F^^mPAH is the 6.2/im PAH flux in units of 
10~ 19 Wcm~ 2 , D the distance in kiloparsecs, and Lir 
in units of Lq. The standard error is 0.3 in log(Lj^). 
Using the above equation, the total infrared luminosity 
of a starburst galaxy can be derived from the strength 
of a single PAH emission feature (here: the 6.2^m PAH) 
to within a factor of two. This correlation is less tight 
than the estimate from the Fis^m and Fso^m continuum 
fluxes (section 4.2), supporting the finding by Peeters et 
al. (2004) that PAHs may not (only) trace recent massive 
star formation. 
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Fig. 14. — The estimated infrared luminosity L^£ H based on 
the aperture corrected 6.2/im PAH flux versus the total infrared 
luminosity as a measure of SFR. The gray plus signs refer to the 
location of the same data points if no aperture correction is applied. 
The dashed line represents a linear fit to the data (in log-log space) , 
and the fitted function is given in the upper left, together with the 
standard error a. The labeled outliers NGC 1365, NGC 4945 were 
excluded from the fit. The filled symbols correspond to starbursts 
with a weak AGN component. 

Combining the information from Figs. 13 and 14 sug- 
gests that the continuum and the PAH emission are, to 
first order, proportional. Similarly, Peeters et al. (2004) 
have found that Lrir is proportional to the PAH lu- 
minosity L 6 ,2fim- If our spectra were "contaminated" 



20 



Brandl et al. 



by a significant amount of continuum emission from an 
AGN or an underlying, older galactic population unre- 
lated to the starburst, the equivalent width would vary 
with distance (equivalent slit width), which is not ob- 
served. Hence, we conclude that both PAH and con- 
tinuum emission originate predominantly from the star- 
burst. 

4.6. PAH Strength and Radiation Field 

It is well known that the equivalent width of PAHs 
is much reduced in AGN-dominated environments (e.g., 
Sturm ct al. (2000), Genzel & Cesarsky (2000), Weed- 
man et al. (2005)). Geballe et al. (1989) have discussed 
the susceptibility of PAHs to destruction by far UV fields, 
and ISO observations (e.g., Cesarsky ct al. (1996), Tran 
(1998)) have shown that more intense far-UV radiation 
fields may lead to gradual destruction of PAHs around 
stellar sources. Recently, Wu et al. (2006) studied a sam- 
ple of low-metallicity blue compact dwarf galaxies from 
1/50Z Q to 1/1.5Z Q with the IRS. They find a strong 
anti-correlation between strength (equivalent width) of 
the PAH features and the product of the [Ne III] /[Nell] 
ratio (as a hardness measure of the radiation field) and 
the UV luminosity density devided by the metallicity. 
A similar trend has been reported by Madden et al. 
(2006) for a small sample of nearby dwarf galaxies. Un- 
fortunately, lower metallicity and harder radiation fields 
seem to go hand-in-hand in these dwarf galaxies, and 
one cannot unambiguously distinguish between possibly 
suppressed PAH formation in low metallicity environ- 
ments and PAH destruction in harder UV fields. Re- 
cently, Beirao et al. (2006) have investigated the strength 
of the 11.3/im PAH feature in the starburst in NGC 5253 
for different radial distances, and found that the equiv- 
alent width of the PAH feature is inversely proportional 
to the intensity of the radiation field, suggesting photo- 
destruction of the aromatic carriers in harsher environ- 
ments. 

Observations of Galactic sources (e.g., Verstraete et al. 
(1996), Vermeij et al. (2002)) have also shown that the 
relative strengths of individual PAH features can depend 
on the degree of ionization of the molecule: C-C stretch- 
ing modes at 6.2/im and 7.7/im are stronger in ionized 
PAHs, while the C-H in-plane bending mode at 8.6/im 
and the C-H out-of-plane bending mode at 11.3/^m are 
stronger by more than a factor of two in neutral PAHs. 
Comparing the ISOCAM spectra of M82, NGC 253, and 
NGC 1808, Forstcr Schreiber et al. (2003) found that, 
while the 5 — 11/im spectrum is nearly invariant, the 
relative PAH intensities exhibit significant variations of 
20% - 100%, which they attributed to the PAH size dis- 
tribution, ionization, dehydrogenation, or the incident 
radiation field. 

In Fig. 15 we compare the ratio between the 11.3/mi 
bending mode and the 7.7/mi stretching mode to the 
hardness of the radiation field, as indicated by the 
[Ne III] /[Nell] ratio from Devost et al. (2006), who pro- 
vide a detailed analysis of the fine structure lines in our 
starburst sample. We correct both PAH fluxes for extinc- 
tion (section 4.4) using the extinction laws by Draine 
(1989) and Lutz (1999) and the values from Table 5. 
Within the systematic uncertainties, represented by the 
scatter of the data points, we find no significant variation 
of the 7.7/zm / 11.3/zm PAH ratio over more than an order 



of magnitude in [Ne III] /[Nell] fluxes. While we cannot 
exclude variations on scales of individual Hn regions, 
our spatially averaged spectra of starburst nuclei do not 
reveal significant variations between the main PAH fea- 
tures. We have also looked at the weaker PAH features 
from Table 4, but the scatter increases with lower signal- 
to-noise and more uncertain baseline definition, and does 
not reveal an obvious trend. However, more studies, in 
particular of the PAH complex around 17/im are planned. 
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Fig. 15. — The ratio of the 7.7/im / 11.3/tm PAH fluxes versus 
the [Ne III] /[Nell] line ratio (from Devost ct al. (2006)). The PAH 
fluxes have been corrected for extinction using the laws from Draine 
(1989) (diamonds) and Lutz (1999) (squares). The gray pluses in- 
dicate the locations of the uncorrected values for reference. The 
filled diamonds correspond to starbursts with weak AGN compo- 
nents. 

Our starburst sample spans a wide range in radiation 
field hardness. Fig. 16 shows the equivalent width of the 
7.7/mi PAH feature versus the fine structure line ratio 
[Ne III] /[Nell], which has been taken from Devost et al. 
(2006). Within the uncertainties the equivalent width of 
the 7.7/im PAH feature remains constant over more than 
an order of magnitude in [Ne III] /[Nell]. We conclude 
that, on large scales of starburst nuclei, which typically 
contain numerous Hn regions, the PAH / continuum ra- 
tio does not significantly depend on the average radiation 
field hardness. 

5. SUMMARY 

We presented and discussed the 5 — 38/im mid-infrared 
spectra of a large sample of 22 starburst galaxies taken 
with the Infrared Spectrograph IRS on board the Spitzer 
Space Telescope. The high signal-to-noise spectra con- 
tain numerous important diagnostics such as PAH emis- 
sion features, silicate bands at 9.8/im and 18/im, and the 
shape of the spectral continuum. The IRS spectral res- 
olution of R « 65 — 130 is perfectly matched to study 
these features. From our sample we constructed an aver- 
age starburst spectrum, which can be used as a starburst 
template. 

The availability of continuous mid-infrared spectra of 
numerous objects within one class over a wide wavelength 
range enables various important studies. Remarkably, 
the spectra show a vast range in starburst SEDs. We 
found a trend that more dust extincted starburst sys- 
tems have a steeper spectral continuum slope longward 
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Fig. 16. — The equivalent width of the 7.7/im PAH feature 
versus the [Ne III] /[Nell] line ratio (from Devost et al. (2006)) 
as an indicator of the hardness of the radiation field. The filled 
diamonds correspond to starbursts with weak AGN components. 
NGC 1365, which has a significant AGN component, is a clear 
outlier. 



of 15/tm. The slope can also be used to discriminate be- 
tween starburst and AGN powered sources, with a transi- 
tion at 0.17 < Fi^m/Fsoum < 0.21. The monochromatic 
continuum fluxes, which represent a more accurate esti- 
mate of the "true" continuum than broad band filters, 
provide a remarkably accurate estimate of the total in- 
frared luminosity via TJf^ = 

(after correcting for slit losses for nearby, extended sys- 
tems). 

Our starburst spectra cover a wide range of silicate ab- 
sorption depths, from essentially no absorption to heav- 
ily obscured systems with an optical depth of Tg.s^m ~ 
5. We present the discovery of crystalline silicates in 
NGC 4945, which shows many similarities with heavily 
extincted ULIRGs. However, unlike the latter, the star- 
bursts in our sample show no signs of water ices or hy- 
drocarbons, suggesting a small amount of self-shielding. 

The observed spectra show significant variations in 
the relative strengths of the individual PAH features at 



6.2/im, 7.7/xm, and 11.3/im. However, these variations 
may be entirely due to extinction and do not necessar- 
ily indicate intrinsic variations of the PAH spectrum. 
We find that the PAH equivalent width is independent 
of the total luminosity L/ij, probably because the PAH 
strength and the underlying continuum scale proportion- 
ally within a "pure" starburst. The luminosity of an 
individual PAH feature, however, scales with Ljr. In 
particular the 6.2/tm feature can be used to approximate 
the total infrared luminosity of the starburst (although 
less accurately than from the 15/um and 30/^m continuum 
fluxes). 

We investigated possible variations of the PAH spec- 
trum as expected, e.g., from varying degrees of PAH ion- 
ization. The 7.7/tm / 11.3/im PAH ratios show no signifi- 
cant systematic variation with the hardness of the radia- 
tion field. Although our sample covers about a factor of 
ten difference in radiation field hardness (as indicated by 
the [Ne III] / [Ne II] ratio) we found no systematic cor- 
relation with the PAH equivalent width. Furthermore, 
we found no systematic differences between "pure" star- 
bursts and galaxies with a weak, non-dominant AGN 
component for most of their spectral properties (except 
for NGC 1365 which shows very weak PAH emission). 

We emphasize that these results are based on spatially 
integrated diagnostics over an entire starburst region. 
Local variations of age, IMF, density or geometry on 
the scales of individual Hil regions or super star clus- 
ters may just average out. However, it is important to 
note that, because of this "averaging out effect" in unre- 
solved sources, starburst nuclei with significantly differ- 
ent global properties may appear as rather similar mem- 
bers of one class of objects. 
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